Excess conductance of a spin-filtering quantum dot 
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The conductance G of a pair of single-channel point contacts in series, one of which is a spin filter, 
increases from 1/2 to 2/3 x e 2 /h with more and more spin-flip scattering. This excess conductance 
was observed in a quantum dot by Zumbiihl et al., and proposed as a measure for the spin relaxation 
time Ti. Here we present a quantum mechanical theory for the effect in a chaotic quantum dot 
(mean level spacing A, dephasing time tq, charging energy e 2 /C), in order to answer the question 
whether Ti can be determined independently of and C. We find that this is possible in a 
time-reversal-symmetry-breaking magnetic field, when the average conductance follows closely the 
formula (G) = (2e 2 /h)(T 1 + /i/A)(4Ti + 3/i/A) -1 . 

PACS numbers: 72.25.Dc,72.25.Rb,73.23.-b,73.63.Kv 
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The study of spin relaxation in the presence of chaotic 
scattering is a challenge for theorists and experimental- 
ists. The common goal is to identify transport properties 
that can be readily measured and that depend as directly 
as possible on the spin relaxation time (T\). One line 
of research is to study how quantum interference effects 
such as weak localization or universal conductance fluc- 
tuations are modified by spin relaxation 1]. A direct 
relation with T± in that context is hindered by the fact 
that dephasing (both of the orbital and of the spin de- 
grees of freedom) also modifies the quantum interference 
effects. Another line of research is to study spin-resolved 
current noise |2| ■ There a direct relation with T\ is possi- 
ble, but the complications involved in the measurement of 
both time- and spin-dependent current fluctuations have 
so far prevented an experimental realization. Ideally, one 
would like to relate T\ to the time averaged current in a 
way which is insensitive to dephasing. It is the purpose 
of this work to present such a relationship. 

Our research was inspired by the proposal of Zumbiihl 
et al. of a new technique to measure spin relaxation times 
in confined systems |3( • These authors reported measure- 
ments of the conductance of an open two-dimensional 
GaAs quantum dot in a parallel magnetic field. One of 
the two point contacts was set to the spin-selective e 2 /h 
conductance plateau. The other point contact was set to 
transmit both spins. In this configuration, the classical 
series conductance of the two point contacts is h x e 2 /h 
if there is no spin relaxation and | x e 2 /h if there is 
strong spin relaxation. What we will show here is that 
the ensemble averaged conductance in a time-reversal- 
symmetry-breaking magnetic field varies between these 
two limits as a rational function of the product of T\ and 
the mean level spacing A — largely independent of the 
presence or absence of dephasing. 

The geometry of the problem is sketched in Fig.Q] We 
discuss its various ingredients. 

Electrons in a two-dimensional electron gas (2DEG) 
enter and leave the quantum dot via two single-channel 
quantum point contacts (QPC). A QPC can operate as 
a spin filter in a magnetic field 0, 01 , as a result of the 
slightly different Fermi wave lengths of spin-up and spin- 
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FIG. 1: Illustration of the model. A current / is passed 
through a quantum dot via two single-channel leads, at a volt- 
age difference V. Spin-flip scattering and decoherence (with 
relaxation times Ti and t</>) are introduced by means of ficti- 
tious voltage probes, separated from the quantum dot by tun- 
nel barriers (dashed lines). The lower (ferromagnetic) voltage 
probe reinjects an electron into the quantum dot with the 
same spin but a random phase (contributing only to t^). The 
upper (normal metal) voltage probe randomizes both spin and 
phase (contributing to both Ti and r^). 



down electrons. The filtering property of a QPC can 
be turned on and off by adjusting its local electrostatic 
potential (via a gate voltage). The polarity of the spin 
filter is fixed by the direction of the magnetic field. The 
conductance becomes sensitive to spin-flip scattering if 
one point contact is a spin filter while the other transmits 
both spin directions. (To be definite, we will take the 
current source as the spin filter, but it does not matter 
which is which in the linear response regime.) 

We assume that the magnetic field is sufficiently weak 
that we may neglect the spin dependence of the Fermi 
wave length inside the quantum dot (where the Fermi 
energy is much greater than in the point contact). The 
effect of the magnetic field on the orbital motion will 
typically break time reversal symmetry (symmetry index 
(3 = 2), if the field is oriented perpendicular to the 2DEG. 
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We contrast this with the case (3 = 1 of preserved time 
reversal symmetry, appropriate for moderately weak par- 
allel fields (until the finite thickness of the 2DEG drives 
p = 1 i— * 2 even for a parallel field [(|). The mean dwell 
time in the quantum dot is assumed to be small compared 
to the spin-orbit scattering time, so that spin-orbit cou- 
pling can be neglected. Landau level quantization inside 
the quantum dot is assumed to be insignificant. The ef- 
fects of a finite charging energy will be assessed at the 
end of the paper. 

Two independent time scales characterize the spin de- 
cay, the time scale T± on which the spin direction is ran- 
domized and the time scale Ti < 2T\ on which the phase 
of the spin-dependent part of the wave function is ran- 
domized @)laEl- I n closed GaAs quantum dots, hyper- 
fine interaction with nuclear spins is the dominant source 
of spin decay for weak magnetic fields, with T2 — /is and 
Ti increasing from ps to ms with increasing magnetic field 
|lCt ll l| . For the transport problem in an open quantum 
dot considered here, the decoherence time of the whole 
wave function, rather than just its spin-dependent part, 
is the relevant quantity. Typically, is dominated by 
dephasing of the orbital degrees of freedom by electron- 
electron interactions. 

Experiments [T^ | on the effect of a finite on spin- 
independent conduction have been analyzed in the past 
using Biittiker's voltage probe model |13lll4lll5| . Exten- 
sions to spin-dependent conduction have been proposed 
more recently [r| El- As described in Ref. |l6j . one 
needs two types of voltage probes to describe spin relax- 
ation and decoherence. One type of voltage probe is con- 
nected to a normal metal reservoir, while the other type 
of voltage probe is connected to a pair of ferromagnetic 
reservoirs (of opposite polarization, parallel to the polar- 
ization of the spin filters in the quantum point contacts). 
For each reservoir, an electron that enters it is reinjected 
into the quantum dot with a random phase. The ferro- 
magnetic reservoirs conserve the spin (contributing only 
to t$), while the normal metal reservoir randomizes the 
spin (contributing both to T\ and r^). 

Each voltage probe is connected to the quantum dot 
by a tunnel barrier. The normal metal voltage probe has 
iV,J = jyl = ]\[ n channels for each spin direction and 
the ferromagnetic voltage probes have Nj — Nj = Nf 
channels. Each barrier has tunnel probability T per chan- 
nel and per spin direction. By taking the limit T — > 0, 
N n ,Nf — > 00 at fixed (dimensionless) tunnel conduc- 
tances 7„ = N n T, 7/ = NfT we ensure that the de- 
cay processes are spatially homogeneous . The decay 
times are 
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T <t> 



(7n + 7/)A 7^A' 



(1) 



with A the mean spacing of spin-degenerate levels and 
70 = 7 n + 7/. These time scales should be compared 
with the spin-dependent mean dwell time T^ well in the 
quantum dot without voltage probes, given by rl wel , = 



The electron reservoirs connected to the quantum dot 
have electrochemical potentials fj, x , with X = s (source), 
X = d (drain), X — n (normal metal voltage probe), 
and X = f (ferromagnetic voltage probe). In the latter 
case we distinguish the two spin polarizations by a super- 
script: (j/fj/j/i. We choose the zero of energy such that 
fid = 0, hence (i s — eV. Both the temperature and the 
applied voltage V are assumed to be small compared to 
A, so that we may neglect the energy dependence of the 
scattering processes. 

The potentials of the voltage probes are determined by 
demanding that no current is drawn from the quantum 

dot PEH, 



= (2N n -T^ n -T^ n )i, n -T^ n eV-T}^y f 

~ T U»P ( 2 ) 
= (N f -Tj^ f )f,}-T^ f eV-T^ n , (3) 

= (N f -Tj^ f )^ f -T^ n . (4) 
The current / through the quantum dot then follows from 



-I = (1 - TU a )eV - Tl^ r , 
e 



(5) 



h/2A, 



h/A. 



Here T X ^ Y and T X ^ Y denote the transmission proba- 
bilities, summed over all channels, from reservoir X to 
reservoir Y with spin up or down. They satisfy the sum 
rules 

E T x^y= E T Y ^ X =N X , (6) 

Y=s,d,n,f Y=s,d,n,f 

with N s — Nd = 1. For later use we define 

^ = 2 - TJL. - T^ d - T^ d - Tj_, -R} = 1- T^ d . 

(7) 

Because of the spatial homogeneity of the coupling of 
the quantum dot to the voltage probes, the transmission 
probabilities for normal and ferromagnetic probes are re- 
lated by ratios of tunnel conductances, 

T x _ f T^ x 7/' 1 h ^ ' 

T X ^ Y -S XY N X {1-T^ B ) = 1X 1Y 

t x ,^ y , - s x . Y -N x , (1 - r- ff ) lx , lY , ' 

if X,Y€ {nj}. (9) 

The effective tunnel probability T° s = (rA)p a differs 
from the bare tunnel probability T because the density 
of states p a in the quantum dot has spin and energy 
dependent fluctuations around the average 1/A. 

With the help of these relations the solution of Eqs. 
for the conductance G — I/V can be written in 
terms of transmission probabilities between source and 
drain, 
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g = \ [i - tJU. - Q(i - rJU. - rj^,)(i - tJU. - r;_j] , (10) 

c = r c T ff r c x ff7n ( 7n + 7/ ) + r(r,[ ff + r^ye* 

r c T ff r* ff7n ( 7n + 7/ )(^i + rci) + r(rT ff + r^TVui ' 



The transmission probabilities between source and 
drain are constructed from two scattering matrices 
and S^, one for spin-up and one for spin-down. The 
spin-up scattering matrix is a 2 x 2 matrix, 



r a 

t r' 



(12) 



such that T] 



T 



„'|2 



T 



■ I' I ' = 1*1: aild 

rj^s = |i'| 2 - The matrix S is symmetric for ft = 1, 
meaning that = 2d_,. s in that case. For ft — 2 the 

two transmission probabilities are not related. Because 
of the voltage probes, S is sub-unitary. The eigenvalues 
t\ , T2 S [0, 1] of the matrix 1 — SfiS^ give the probability 
to enter one of the voltage probes. 

The statistics of the matrix Sr in an ensemble of 
chaotic quantum dots was calculated in Ref. [l5] using 
the methods of random-matrix theory. It is given in 
terms of the polar decomposition 



S 







(13) 



with unitary matrices vl = u T if ft = 1 and v! inde- 
pendent of u if ft = 2. These matrices are uniformly dis- 
tributed in the unitary group. The distribution Pp(r\, 72) 
is the Laguerre ensemble for 7^, <C 1 and a more compli- 
cated (but known) function for larger 7^. 

In addition to we also need . This is a single com- 
plex number, such that T^ d = {S^ 2 . It is constructed 
from the coefficients r, t, t' in Eq. (|13fl by reflecting spin- 
down from the source contact, 



S l = r' 



Ht' 



1 



(14) 



(The phase shift a need not be specified because it drops 
out upon averaging over u and u'.) Using Eq. i|14|) the 
statistics of S 1 ^ follows from the statistics of S' . 

To complete the random-matrix theory, we need to 
know the statistics of the density of states p a of the open 
quantum dot, which determines the effective tunnel prob- 
abilities. For weak decoherence we have the relation Oil 



Tr(l-S CT S CT t) = ^ 70 A + 0( 7 2). 



(15) 



Since the left-hand-side of Eq. (|15l) equals Vf by defini- 
tion |(7|), we have 



In the opposite limit 7^ > 1 the fluctuations in the den- 
sity of states can be neglected, so that 



r e vr = p ct a = 1, if 14> » 1. 



(17) 



These two limits are sufficient for the purpose of compar- 
ing coherent and incoherent regimes. 

We calculate the average conductance (G) separately 
in the regime 7 / ^> 1 of strong orbital dephasing and 
the regime 7/ <C 1 of weak orbital dephasing. For strong 
dephasing we have T° e — ► T, IZ^ — » 2, V} — ► 1, t\,tz — > 
1, hence 



2e 2 1 + 7n 

(G) = -r- . ~ o , lf 7/ > 1- 

h 4 + 3 7 „ 



(18) 



This incoherent regime is insensitive to the presence or 
absence of time-reversal symmetry. By writing Eq. 1|18|) 
as (G) = (e 2 // l )[2(l-p) + |p]- 1 withp = 7n /(l+ 7n ), we 
can understand it as a classical series resistance, weighted 
by the probability p of a spin- flip scattering event. 

Turning now to the phase coherent regime, we find to 
linear order in 7/ and 7 „ the expansions 

-(C) - [ 3 + °- 147 " + if & = lj (Vt\ 

e 2[ ) ~\ i + 0.10 7n , if ft = 2. [iJ) 

Note the absence of a term linear in 7/ for ft = 2. The 
difference between the zeroth order terms 1/3 and 1/2 
in the presence and absence of time-reversal symmetry is 
known as weak localization or coherent backscattering. 

In the absence of any orbital dephasing, 7/ = 0, we 
obtain the results plotted in Fig. [21 (dashed curves). 
Comparison with the incoherent result (|18|) (solid curve) 
shows that the presence or absence of orbital dephasing 
does not change (G) by more than a few % if ft = 2 
(no time-reversal symmetry). For ft = 1, in contrast, the 
dependence on T\ is entirely different with and without 
orbital dephasing. 

So far we have not included the effects of a finite charg- 
ing energy e 2 /G (with C the capacitance of the quan- 
tum dot). These results therefore apply to the regime 
e 2 /G -C A. In the opposite, more realistic, regime 
e 2 /G 3> A the charging energy introduces a weight fac- 
tor equal to the density of states in the ensemble averages 
[l^| . This weight factor converts the grand-canonical av- 
erage (• • •) (considered so far) into a canonical average: 



r: ff /r = p °a = n°/^, if 70 « 1. 



(16) 



/canonical 



\M---(p ] +p 1 ))- 



(20) 
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FIG. 2: (color online) Dependence of the average conduc- 
tance (G) on the spin relaxation time Ti, normalized by the 
mean level spacing A. The solid curve (red) is the incoherent 
result (1181 . valid for strong orbital dephasing. The dashed 
and dotted curves are the results of random-matrix theory 
for weak orbital dephasing, in the two cases of broken (/3 = 2, 
black) and unbroken ((3—1, blue) time reversal symmetry. 
The dashed curves are grand-canonical averages (e 2 /C -C A) 
and the dotted curves are canonical averages (e 2 /C 3> A). 
The black and red curves lie close together, demonstrating 
that Ti can be determined accurately from (G) for (3 = 2. 
The large difference between the blue and red curves prevents 
this for 13 = 1. 



The incoherent result l|18|) is the same in the canoni- 
cal and grand-canonical ensembles, because the fluctu- 
ations in the density of states are suppressed by deco- 
herence. In order to assess the importance of density-of- 
states fluctuations in the coherent regime, we approxi- 
mate (A/2)(pT + pl) ~ (ftT +TZ l )/(TZ^ +K l ). This for- 
mula interpolates smoothly between the two exact limits 
(|16fl and (|17fl of weak and strong decoherence. As shown 
in Fig. (dotted curves), the effect on the average con- 
ductance remains relatively small. 



In conclusion, we have calculated the dependence on 
the spin relaxation time Tj of the average conductance 
(G) of a quantum dot with a spin-filtering quantum point 
contact. In the incoherent regime there is a simple one- 
to-one relationship (|18fl between (G) and T\. The pres- 
ence or absence of orbital dephasing was found to be 
insignificant for f3 = 2, so that the value of T\ can be 
extracted from (G) with good accuracy — without re- 
quiring knowledge of coherence time or charging energy. 
For P = 1, in contrast, the interplay with the weak local- 
ization effect obscures the effect of spin relaxation. 
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